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ABSTRACT 

Aims. We extend our previous work on the effects of the uncertainties on the main input physics for the evolution of low-mass stars. 
We analyse the dependence of the cumulative physical uncertainty affecting stellar tracks on the chemical composition. 
Methods. We calculated more than 6000 stellar tracks and isochrones, with metallicity ranging from Z = 0.0001 to 0.02, by changing 
the following physical inputs within their current range of uncertainty: 'H(p,ve + ) 2 H, 14 N(p, y) 15 and triple-o- reaction rates, radiative 
and conductive opacities, neutrino energy losses, and microscopic diffusion velocities. The analysis was performed using a latin 
hypercube sampling design. We examine in a statistical way - for different metallicities - the dependence on the variation of the 
physical inputs of the turn-off (TO) luminosity, the central hydrogen exhaustion time (f H ), the luminosity and the helium core mass at 
the red-giant branch (RGB) tip, and the zero age horizontal branch (ZAHB) luminosity in the RR Lyrae region. 
Results. For the stellar tracks, an increase in the metallicity from Z = 0.0001 to Z = 0.02 produces a cumulative physical uncertainty 
error variation in TO luminosity from 0.028 dex to 0.017 dex, while the global uncertainty on tn increases from 0.42 Gyr to 1.08 
Gyr. For the RGB tip, the cumulative uncertainty on the luminosity is almost constant at 0.03 dex, whereas the one the helium core 
mass decreases from 0.0055 M Q to 0.0035 M Q . The dependence of the ZAHB luminosity error is not monotonic with Z, and it varies 
from a minimum of 0.036 dex at Z = 0.0005 to a maximum of 0.047 dex at Z = 0.0001. Regarding stellar isochrones of 12 Gyr, the 
cumulative physical uncertainty on the predicted TO luminosity and mass increases respectively from 0.012 dex to 0.014 dex and from 
0.0136 Mq to 0.0186 M Q . Consequently, fromZ = 0.0001 toZ = 0.02 forages typical of galactic globular clusters, the uncertainty on 
the age inferred from the TO luminosity increases from 325 Myr to 415 Myr. 

Key words, methods: statistical - stars: evolution - stars: horizontal-branch - stars: interiors - stars: low-mass - stars: Hertzsprung- 
Russell and C-M diagrams 
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\£> 1. Introduction helium-core mass M^ e at the red giant branch (RGB) tip, and 

t^ . i 1 the ZAHB luminosity Lhb in the RR Lyrae region at log T e s = 

• The present paper extends the analysis presented in |Valleetal. | 3 . 83 . A s turn-off luminosity we adopted the luminosity of a point 

g (|20L3b (hereafter PI) on the impact of the physical uncertainties brighter and 100 K lower than the turn- off (hereafter BTO), a 

^ affecting stellar evolution models to a wide range of chemical technique similar to the one proposed bv lChaboveret"aII(lL996l) 

,_, compositions covering the typical values for stellar populations for isocnrones i n the (B-V, V) plane and already adopted in PI. 
.. in the Milky Way and in the near galaxies. As in the previous Table [T]lists the physical input that was allowed to vary and 

m £ work, this paper is focused on the evolutionary characteristics of the assumed uncertainty on each. We refer to PI for a detailed 

k> ancient stellar populations, and thus the analysis is restricted to discussion underlying their choice and the determination of the 

S_j stellar models of low-mass stars from the main sequence (MS) uncertainty ranges 
03 to the zero age horizontal branch (ZAHB). 

In PI we determined, for a fixed chemical composition 
(metallicity Z = 0.006, and helium abundance Y - 0.26), the 



2. Description of the technique 



range of variation in some important evolutionary features due to We adopt a reference model of M = 0.90 M to explore the cu- 
the changes in the input physics within their uncertainty, point- mulative effect of the uncertainty on the main input physics, for 
ing out the effect of the different physical inputs at different s i x different metallicity values, i.e. Z = 0.0001, 0.0005, 0.001, 
stages of stellar evolution. 0.006, 0.01, and 0.02. 

The aim of this paper is to analyse how the results found Determination of the initial helium abundances was done us- 

in PI depend on the assumed chemical composition. Thus, we ing a linear helium-to-metal enrichment law: Y - Y p + ^|Z, with 
computed models covering a wide metallicity range, i.e. from Z cosmological 4 He abundance Y D = 2485 dCvburt et al.ll2004t 
= 0.0001 to 0.02. We focus on the turn-off luminosity L BT o, the Isteigm an 2006; iPeimbert et ai1l2007allbl) . We adopted AY/AZ = 
central hydrogen exhaustion time t H , the luminosity L tip and the 2 , as the typical value for this quantity, whi ch is still affected 

by several important sources of uncerta inty (Pagel & Portinari 
119981: iJimenez et all l2003t iFlvnnl 1200 4: Genn aroetall EoiO). 
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Table 1. Physical inputs perturbed in the calculations and their assumed 
uncertainty. The abbreviations used in the following tables are defined 
in parentheses. 



description 


parameter 


uncertainty 


'H(p,ve + ) 2 H reaction rate (pp) 


Pi 


3% 


14 N(p,y) ls O reaction rate ( 14 N) 


Pi 


10% 


radiative opacity (k r ) 


P3 


5% 


microscopic diffusion velocities (vd) 


IH 


15% 


triple-a reaction rate (3a) 


P5 


20% 


neutrino emission rate (v) 


P6 


4% 


conductive opacity (k c ) 


Pi 


5% 



Notes. The table is taken from PI and is reported for the reader's con- 
venience. 



Thus, we computed models by adopting the following couples of 
initial metallicity and helium abundance (Z, Y): (0.0001, 0.249), 
(0.0005, 0.25), (0.001, 0.25), (0.006, 0.26), (0.01, 0.268), and 
(0.02, 0.288). 

The adopted stellar evolutionary code, FRANEC, is the same 
as used in PI and for the construct ion of the Pisa Stellar Evolu - 
tion Data Basq^for low-mass stars (DeH'Omodarme et al. 2012). 
We refer to these papers and to Valle et al. (2009) for a detailed 
description of the input on the stellar evolutionary code and of 
the ZAHB construction technique. 

We briefly recall the technique followed in PI, i.e. a system- 
atic variation in the input physics on a fixed grid within their 
current uncertainty. For every physical input listed in Table Q] 
we adopted a three-value multiplier p-, with value 1 .00 for the 
unperturbed case and values 1 .00 + A/?, for perturbed cases (A/?, 
is the uncertainty listed in Table[TJ- For each stellar track calcu- 
lation, a set of multiplier values (i.e. p\,... ,pi) is fixed. Then, 
to cover the whole parameter space and explore all the possible 
interactions among input physics, calculations of stellar tracks 
must be performed for a full crossing, i.e. each parameter value 
Pi crossed with all the values of the other parameters pj, with 
;' + i. 

However, as shown in PI, the interactions among the physical 
inputs are negligible for such small perturbations. This makes 
possible to reduce the computational burden by avoiding compu- 
tation of the whole set of 3 7 = 2187 tracks with the same mass, 
chemical composition, and a m \ for each metallicity. We select 
randomly - for each Z - a subset of n — 162 models to be com- 
puted using a latin hypercube sampling design: it is an extension 
of the latin square to higher dimensions and it has optimal prop- 
erty in reducin g the varian ce of the estimators obtained from the 
linear models (Stein 1987). The sample size was chosen to bal- 
ance the expected run time with good precision in the estimation 
of the effects of the input physics. In particular, we evaluated 
the impact of a sample size reduction on the significance of the 
estimated coefficients reported in PI. For this purpose, we con- 
sidered that the estimated errors of the effects (see Sect. [3} of 
the input physics - hence the confidence intervals width of the 
effects - scale approximately as the square root of the ratio of 
the original sample size (i.e. 2187) and the subsample size n. 

Test simulations, performed by extracting different samples 
of different size n and analysing the influence of the input 
physics in each sample, showed that the sample size n — 162 is 
sufficient to obtain robust estimates of the input physics impacts. 
The random selection was performed using the R library lhs 
(Carnell 2012), resulting in a total sample of 972 stellar tracks. 



3. Statistical analysis of physical uncertainty 

Relying on such a large set of perturbed stellar models, we statis- 
tically analysed the effects of varying the physical inputs listed in 
TableQ]on the selected evolutionary features. For each Z value, 
we built linear regression models by extracting the values of the 
dependent variable under study (i.e. Lbto, ?h, L t ; p , M^ e , and 
Lhb) from the stellar tracks and regressing them against the co- 
variates, i.e. the values of the parameter /?,. Thus we stratified 
the sample according to Z values and we treated strata separately. 

An alternative approach would require a global fit of the data, 
modelling the dependence on Z and Y by inserting of interac- 
tion terms into the statistical model; however, direct computa- 
tions showed that the models of these dependencies are quite 
cumbersome and require a function with several powers of log- 
metallicity. We therefore prefer not to rely on these complex 
models and to put our trust in the simpler ones obtained by strat- 
ification. 

The statistical models include only the predictors but no in- 
teractions among them, since they can be safely neglected, as 
we showed in PI. For each studied evolutionary feature, we in- 
cluded as covariates in the model only the physical input that 
can have an actual influence: for Lbto an d f H the first four pa- 
rameters of Table Q] and all the parameters for the quantities re- 
lated to later evolutionary stages. The models were fitted to the 
data with a least-squares method using the software R 2.15.1 
(IR Development Core Teamll2012|) . 

For each Z value, for log Lbto and f H the regression models 
are 



log Lbto, t H = yS + £jPi Pi 



(1) 



where /?o, ■ ■ • ,Pa were the regression coefficients to be estimated 
by the fit and p\,...,pi the perturbation multipliers defined in 
Table[T] 

In the cases of logL t ; p , log Lhb, and M^ e , the linear models 
for each Z value are 



logLtip, log L HB , M c He = Po + ^ A Pi 



(2) 



1 http://astro.df.unipi.it/stellar-models/ 



where we added the dependence on the triple-a reaction rate, the 
neutrino emission rate, and the conductive opacity k c . 

The cumulative effect of the physical input perturbation ob- 
tained from the statistical models for each evolutionary feature 
is listed in Table [2] There we report in the second column the 
reference value of the studied quantity obtained from the model 
with unperturbed physics, for each Z value, in the third column 
the cumulative statistical impact of the various physical inputs 
defined as £,■ A/?,- x |/?,|. 

>From Table|2]some general trends are apparent. In the case 
of BTO log-luminosity, the total impact decreases with metallic- 
ity from Z = 0.0001 to 0.02 of about 40%. For central hydrogen 
exhaustion time, we note that the total impact grows with the 
metallicity with a rate slightly higher than the one of the increase 
in the reference time. In the case of RGB tip, we observe that 
the total impact on the luminosity is nearly constant at about 0.03 
dex, whereas on the helium core mass it decreases from 0.0055 
M Q at Z = 0.0001 to 0.0035 M at Z = 0.02. In the case of 
ZAHB luminosity, the effect of the chemical composition on the 
cumulative physical uncertainty is not monotonic. In this case 
the total impact shows an initially high value at Z = 0.0001, a 
sudden drop at Z = 0.0005, and a subsequent increase with the 
metallicity. 



Article number, page 2 ofllOl 



Valle, G. et al.: Cumulative physical uncertainty in modern stellar models (RN) 



o 

o - 














O " 



















O 


O 


O 


° 


° 


o - 

O _ 


X 
A 


A 


A 


X 
A 


X 

A 


A 


p - 

7 . 














o _ 








+ 


+ 


+ 


o - 
l _ 




+ 


+ 










+ 












M - 

O - 















o 



<? ° 



o 'H (p, ve + ) 2 H 










A 14 N (p, y) 15 
+ k t 








+ 


x Diffusion 












+ 






+ + + 








■* §--| 


o 

-A- 

X 


"A" 


o 
""A" 



-4.0 



-3.5 



-3.0 -2.5 
logZ 



-2.0 



-1,5 



-4.0 



-3.5 



-3.0 -2.5 
logZ 



-2.0 



-1.5 





, 1 , , 












O 3a 








V Neutrinos 


V 
A 


V 
A 


V 
A 


B ^c 


i « a 


X 

"0 


X 

Cr " 


X 

- o ' 


5--2---S--- 


H 


IS 


H 


B h H 





o 


O 


5 o o 


+ 










+ 


+ 


t + + 



V 


¥ 


¥ 








+ 


X 


X 


X 


¥ 

X 


? 


"0 


o 





"0" 


" O" 


a 


A 


A 


A 

H 


A 
B 


A 


A 

H 


H 


H 






















o 




o 


o 








o 













s 



V 










A 






.&____ 


8 

o- 


-X- 


_S.__ 


i 


£ 


§ 


B 


X 


X 


2 


o 


© 














O 





O 




+ 


+ 


+ 


+ 


+ 


+ 













-1.5 



Fig. 1. 

values. 



log Z log Z log Z 

Physical impacts - as in Tables lA.3IA.7l - of the increment Ap, of each multiplier p, for the considered physical inputs at different 



loaZ 



In Fig. Q] we display the results of the fits for the different 
evolutionary quantities. This figure shows the physical impact of 
the increment Ap, on each multiplier /?,■ in dependence on log Z. 
The full results of the fits are reported in Tables lA31IA.7l (avail- 
able online). These tables list, for different values of Z, the least 
squares estimations of the regression coefficients /?,, their stan- 
dard error, the physical impact of the increment of Ap, on each 
multiplier /?,, and the absolute value of the relative importance 
of these impacts with respect to the most important one. 

In the case of Z = 0.006 the values of the estimated coeffi- 
cients are different from the ones reported in PI, since the present 
estimates are obtained from a subsample of the original dataset. 
As expected, in almost all cases the present estimates and the 
ones of PI agree within their errors. 

The figure and the tables show that the effect of perturbing 
the main physical input within their current uncertainty on the 
studied stellar evolutionary quantities changes in a rather com- 
plex way when varying the chemical composition. In particu- 
lar, the relative importance of the various input physics in con- 
tributing to the global uncertainty changes with metallicity. As 
a consequence, one should be careful to extrapolate the results 
we obtained in PI for stellar models of Z = 0.006 to significantly 
different metallicities. 

In the case of BTO luminosity, from Fig. [T]and from the last 
column of Table IA.31 we note that the effect of the variation in 
radiative opacity is dominant for each value of Z, while the effect 
of the other physical inputs change mildly with metallicity. The 
second most important input is the one due to the 14 N(p,y) 15 
reaction rate at each Z. The impact of the perturbations scales 
almost linearly with log Z. 



For central hydrogen exhaustion time in Table IA.4I it is ap- 
parent that the importance of radiative opacity is dominant for 
each Z and that the contribution of microscopic diffusion and of 
'H(p,ve + ) 2 H reaction rate increase with metallicity. At Z = 0.02 
they jointly account for about 30% of the impact of k r . Figure Q] 
shows that the impact of the radiative opacity variation grows 
more than linearly with log Z. 

In the case of RGB tip luminosity, radiative opacity is the 
most important uncertainty source, but at lower Z its importance 
is close to the one of triple-a (87%) and neutrino cooling (67%); 
in contrast, for Z = 0.02 the cumulative impact of all the other 
sources of uncertainty is about 80% of the one of k r . Figure [T] 
shows that the effect is due to the growth in the importance of 
the radiative opacities uncertainty with log Z, with a concurrent 
decrease in the importance of the other sources of uncertainty. 

For the helium core mass at RGB tip, the triple-a reaction 
rate is the most important uncertainty source. A distinctive fea- 
ture is the decrease in the importance of neutrino cooling with 
increasing Z; while it accounts for 72% of the impact of the 
triple-a at Z = 0.0001, this percentage decreases to 33% at Z 
= 0.02. We note also from Fig. Q]a non-monotonic trend with 
log Z for the radiative opacities effect. 

For ZAHB luminosity, the radiative opacity is the main 
source of uncertainty at each metallicity. We observe from Fig.Q] 
an increase in the importance of triple-a moving toward high Z. 
A distinctive feature is the drop in the absolute value of the im- 
portance of radiative opacity from Z = 0.0001 to Z = 0.0005, 
followed by a mild increase with logZ. The trend in the cu- 
mulative uncertainty is further complicated since the impact of 
14 N(p,y) 15 0, which is negative at Z = 0.0001, becomes positive 
at higher Z. Moreover, each physical input variation affects both 
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Table 2. Cumulative impact of the uncertainty on the physical inputs on 
the considered evolutionary features. 



z 


Reference 


Impact 


log Lbto (dex) 


0.0001 


0.749 


0.028 


0.0005 


0.648 


0.024 


0.0010 


0.584 


0.023 


0.0060 


0.354 


0.020 


0.0100 


0.270 


0.019 


0.0200 


0.151 


0.017 


ta (Gyr) 


0.0001 


7.21 


0.42 


0.0005 


7.38 


0.45 


0.0010 


7.74 


0.48 


0.0060 


10.54 


0.71 


0.0100 


12.32 


0.84 


0.0200 


15.23 


1.08 


log L t ip (dex) 


0.0001 


3.257 


0.032 


0.0005 


3.315 


0.031 


0.0010 


3.342 


0.031 


0.0060 


3.406 


0.030 


0.0100 


3.419 


0.029 


0.0200 


3.431 


0.029 


M c He (M ) 


0.0001 


0.5023 


0.0055 


0.0005 


0.4936 


0.0053 


0.0010 


0.4909 


0.0048 


0.0060 


0.4836 


0.0041 


0.0100 


0.4807 


0.0037 


0.0200 


0.4750 


0.0035 


log L HB 


(dex) at log r eff = 3.83 


0.0001 


1.792 


0.047 


0.0005 


1.701 


0.036 


0.0010 


1.671 


0.037 


0.0060 


1.565 


0.041 


0.0100 


1.522 


0.043 


0.0200 


1.465 


0.044 



Notes. First column: the metallicity of the stellar model; second col- 
umn: the reference value for the model with unperturbed physical input; 
third column: the sum of the impacts of the analytical fit. 



the He-core mass (and thus indirectly Lhb) and Lhb directly. 
These effects are the source of the non-monotonic trend on the 
cumulative uncertainty shown in Table [2] The precise under- 
standing of the effect of each physical input is made even more 
difficult by the ZAHB luminosity being evaluated at fixed r e ff 
for all Z. This implies that at low Z the region is populated with 
higher masses with respect to the ones at high Z. 



4. Physical uncertainty on stellar isochrones 

The previous statistical analysis of the physical uncertainties af- 
fecting stellar tracks can be extended to isochrones. We focus 
here on isochrones with ages in the range 8-14 Gyr, which is 
suitable for galactic globular clusters. As in PI, we are inter- 
ested in studying the variation near the turn-off region; so we 
performed calculations by varying only the physical input that 
can influence the evolution until this phase, i.e. p\, . . . ,p^. For 
each Z value, we selected stellar models from the full grid by 
varying p\,...,p\ (i.e. 3 4 = 81 cases). Each metallicity grid 



contains several stellar models with different masses, chosen to 
accurately reconstruct - in the desired age range - the BTO re- 
gion. The computed stellar models are listed in Table IA.8I A 
total of 5832 models and 3402 isochrones were calculated. 

Following the technique outlined in PI, we pool together the 
results obtained from isochrones of different ages. For each Z we 
first remove the trend due to age by adapting an ANOVA model 
to the turn-off luminosity log Lg° and mass M^° , with age as 
a categorical predictor. This method increases the power of the 
analysis due to the larger sample examined after pooling. The 
residuals of these models are then regressed against p\, ... ,p\. 
Further details on the technique and discussion of the underlying 
hypotheses can be found in PI. 

The results about the total uncertainty on log Lg° Q and M^° Q , 
due to the simultaneous variation in all the considered input 
physics, are shown in Table [9] As in Sec.|3]the total uncertainty 
is obtained by adding the effect of each covariate. The total un- 
certainty on the reconstructed log Lg° and on the Mg° shows a 
moderate increase with Z. In the case of log Lg° the difference 
between maximum and minimum uncertainty is of the order of 
20%, whereas for M™ Q it is about 37%. 

Table 9. Cumulative impact on the uncertainty on isochrone BTO lu- 
minosity (dex) and on the mass at isochrone BTO (M Q ) for different 
metallicity values. 



Impact 


Z 


log L^ (dex) M™ n (M Q ) 


0.0001 


0.0118 


0.0136 


0.0005 


0.0121 


0.0141 


0.0010 


0.0122 


0.0144 


0.0060 


0.0129 


0.0164 


0.0100 


0.0138 


0.0174 


0.0200 


0.0143 


0.0186 



Notes. The results are obtained by pooling together isochrones of ages 
in the range [8 - 14], see text for details. 



Figure|2]and Tables fTOlfTTI display the detailed results of the 
fits for logLg° Q and M l *° Q . The figure displays the physical 
impact of the increment of A/?, of each multiplier p-, in depen- 
dence on log Z. The tables report, for different Z values, the least 
squares estimations of the regression coefficients /?,, their stan- 
dard error, the physical impacts of the increment of Ap,- of each 
multiplier /?,, and the absolute value of the relative importance 
of these impact with respect to the most relevant one. 

For Mg° (Table fTTTi the relative contribution of the physical 
sources of uncertainty is almost the same for all the metallicity 
values; the only exception is the role of the p-p chain whose rela- 
tive contribution to the cumulative uncertainty slightly increases 
with metallicity. The radiative opacity dominantes over all the 
other physical input. 

In the case of log Lg° (Table IT0T > the dependence on the 
chemical composition shows much more complex behaviour. At 
Z = 0.0001 the main contribution to the uncertainty comes from 
the microscopic diffusion velocities. At Z = 0.001 the effect of 
microscopic diffusion velocities, 14 N(p,y) ls O reaction rate, and 
radiative opacity are similar; whereas at larger Z the effect of k r 
becomes dominant. As shown in Table [9] these variation nearly 
compensate one another in the whole analysed range of metal- 
licity, so that the cumulative uncertainty on log Lg° is nearly 
constant and ranges from 0.012 dex to 0.014 dex. 

In Table [T2l we present the results of the analysis on the un- 
certainty affecting the age inferred from BTO luminosity. We 
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Fig. 2. Physical impacts on isochrone BTO luminosity and mass at isochrone BTO due to the increment Ap, of each multiplier p t for the 
considered physical inputs at different values of log Z. 



constructed - for each Z - linear models linking log L™ Q and 
the logarithm of the age. We ignore the dependence on the phys- 
ical inputs, confounding their effect with the residual standard 
error. The models are then used to perform a reverse inference 
on the value of the age, given L l ^, Q . For each metallicity set, we 
fix the reference value of log L"° Q to the one reached at 12 Gyr. 
More details about the technique are given in Appendix B of PI. 
The results collected in Table fT2lshow that the uncertainty in the 
inferred age increases by almost 30% from Z = 0.0001 to Z = 
0.02, ranging from 325 Myr to 415 Myr. 



5. Conclusions 



In this paper we extended the work presented in IValle et all 
(l2013h that analyses the cumulative uncertainty due to physical 
inputs in stellar models of low-mass stars. The aim was to anal- 
yse the influence on the results of a variation in the assumed 
chemical composition. Thus we calculated stellar models for a 
wide range of metallicity from Z = 0.0001 to Z = 0.02. 

For the stellar tracks, the cumulative physical uncertainty 
in the analysed quantities changes with the chemical compo- 
sition. For turn-off luminosity the cumulative uncertainty de- 
creases from 0.028 dex at Z = 0.0001 to 0.017 dex at Z = 0.02. 
The global uncertainty on the central hydrogen exhaustion time 
increases with Z from 0.42 Gyr to 1.08 Gyr, with a rate that is 
slightly higher than the increase in the reference time of hydro- 
gen exhaustion at any given Z. In the case of the RGB tip, the cu- 
mulative uncertainty on the luminosity is almost constant at 0.03 
dex, while on the helium core mass it decreases from 0.0055 M Q 
at Z = 0.0001 to 0.0035 M at Z = 0.02. The dependence on the 
metallicity of the error in the predicted ZAHB luminosity is not 
monotonic. The total uncertainty shows a value of 0.047 dex at 
Z = 0.0001, a sudden drop to a value of 0.036 dex at Z = 0.0005, 
and a following increase with Z up to a value of 0.044 dex at Z 
= 0.02. 

We confirm the results presented in PI; i.e., for all the anal- 
ysed stellar tracks features (except the He-core mass at the red 
giant branch tip), in the full adopted range of chemical composi- 
tions, the effects of the uncertainty on the radiative opacity tables 
is dominant. The uncertainty on the helium core mass at RGB 
tip is instead mainly due to the variation in the triple-a reaction 
rate. In conclusion, an increase in the precision of the radiative 
opacity tables is mandatory to further improve stellar evolution 
theoretical predictions. 



For the stellar isochrones of 12 Gyr, the cumulative uncer- 
tainty on the predicted turn-off luminosity and mass show a mod- 
erate increase with Z. From Z = 0.0001 to 0.02, the global phys- 
ical uncertainty increases of the order of 20% for log L l ^, Q and 
37% for Mg° Q . As a consequence, for ages in the range 8-14 
Gyr, the uncertainty on the age inferred from the turn-off lumi- 
nosity increases by about 30% from Z = 0.0001 to Z = 0.02, 
ranging from 325 Myr to 415 Myr. 
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Appendix A: On the EOS influence 

As discussed in PI, the uncertainty on the equation-of-state can- 
not be treated with the approach applied to the other physical 
input. To provide a rough estimate of the impact of the present 
uncertainty on EOS, we computed the reference tracks - i.e. the 
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Table 10. Effect of the various physical inputs on isochrones BTO lu- 
minosity (dex), for different Z values. The column legend is the same 
as in TablelA.3l 



Table 11. Effect of the various physical inputs on the masses at 
isochrones BTO (M Q ), for different Z values. The column legend is 
the same as in TablelA.3l 





Estimate 


Std. Error 


Impact 


Importance 




Estimate 


Std. Error 


Impact 


Importance 






Z = 0.0001 










Z = 0.0001 






Po 


-0.0591 


0.0046 






Po 


-0.2388 


0.0019 






Pi (PP) 


0.0775 


0.0038 


0.0023 


0.573 


Pi (PP) 


0.0165 


0.0016 


0.0005 


0.042 


Pi ( 14 N) 


-0.0332 


0.0011 


-0.0033 


0.817 


Pi ( 14 N) 


-0.0019 


0.0005 


-0.0002 


0.016 


Pi (k r ) 


0.0418 


0.0023 


0.0021 


0.515 


Pi (kr) 


0.2328 


0.0009 


0.0116 


1.000 


Pa (Vd) 


-0.0271 


0.0008 
Z = 0.0005 


-0.0041 


1.000 


Pa (v d ) 


-0.0086 


0.0003 
Z = 0.0005 


-0.0013 


0.110 


Po 


-0.0919 


0.0041 






Po 


-0.2480 


0.0015 






Pi (PP) 


0.0893 


0.0034 


0.0027 


0.795 


Pi (PP) 


0.0207 


0.0012 


0.0006 


0.052 


Pi ( I4 N) 


-0.0319 


0.0010 


-0.0032 


0.945 


Pi ( 14 N) 


-0.0023 


0.0004 


-0.0002 


0.019 


Pi (k T ) 


0.0569 


0.0020 


0.0028 


0.844 


Pi (kr) 


0.2382 


0.0007 


0.0119 


1.000 


Pa (Vd) 


-0.0225 


0.0007 
Z = 0.001 


-0.0034 


1.000 


Pa (v d ) 


-0.0086 


0.0002 
Z = 0.001 


-0.0013 


0.108 


Po 


-0.0983 


0.0038 






Po 


-0.2558 


0.0015 






Pi (PP) 


0.0896 


0.0031 


0.0027 


0.826 


Pi (PP) 


0.0238 


0.0012 


0.0007 


0.059 


Pi ( I4 N) 


-0.0325 


0.0009 


-0.0033 


1.000 


Pi ( 14 N) 


-0.0025 


0.0004 


-0.0002 


0.020 


Pi (k T ) 


0.0622 


0.0018 


0.0031 


0.957 


Pi (h) 


0.2432 


0.0007 


0.0122 


1.000 


Pa (Vd) 


-0.0210 


0.0006 
Z = 0.006 


-0.0031 


0.968 


Pa (v d ) 


-0.0087 


0.0002 
Z = 0.006 


-0.0013 


0.107 


Po 


-0.1400 


0.0027 






Po 


-0.2998 


0.0021 






Pi (PP) 


0.1038 


0.0022 


0.0031 


0.747 


Pi (PP) 


0.0413 


0.0017 


0.0012 


0.092 


Pi ( 14 N) 


-0.0295 


0.0007 


-0.0029 


0.707 


Pi ( 14 N) 


-0.0024 


0.0005 


-0.0002 


0.018 


Pi (kr) 


0.0834 


0.0013 


0.0042 


1.000 


Pi (k r ) 


0.2702 


0.0010 


0.0135 


1.000 


Pa (v d ) 


-0.0177 


0.0004 
Z=0.01 


-0.0026 


0.635 


Pa (v d ) 


-0.0093 


0.0003 
Z = 0.01 


-0.0014 


0.103 


Po 


-0.1545 


0.0031 






Po 


-0.3200 


0.0015 






Pi (PP) 


0.1112 


0.0026 


0.0033 


0.725 


Pi (PP) 


0.0493 


0.0012 


0.0015 


0.104 


£2 ( 14 N) 


-0.0293 


0.0008 


-0.0029 


0.637 


Pi ( 14 N) 


-0.0023 


0.0004 


-0.0002 


0.016 


Pi (kr) 


0.0920 


0.0015 


0.0046 


1.000 


Pi (k r ) 


0.2833 


0.0007 


0.0142 


1.000 


Pa (Vd) 


-0.0195 


0.0005 
Z = 0.02 


-0.0029 


0.635 


Pa (v d ) 


-0.0103 


0.0002 
Z = 0.02 


-0.0015 


0.109 


Po 


-0.1738 


0.0031 






Po 


-0.3438 


0.0018 






Pi (PP) 


0.1164 


0.0025 


0.0035 


0.673 


Pi (PP) 


0.0576 


0.0014 


0.0017 


0.115 


jS 2 ( 14 N) 


-0.0273 


0.0008 


-0.0027 


0.526 


Pi ( 14 N) 


-0.0016 


0.0004 


-0.0002 


0.011 


/?3 (fe) 


0.1038 


0.0015 


0.0052 


1.000 


Pi (kr) 


0.2994 


0.0009 


0.0150 


1.000 


A (Vd) 


-0.0191 


0.0005 


-0.0029 


0.553 


Pa (v d ) 


-0.0115 


0.0003 


-0.0017 


0.115 



ones with unperturbed physics - with two different and widely 
adopted ch oices for equation-of-state: OPAL EOS and FreeEOS 
<IIrwinll2004l) . 

Table lA. ll reports - for each evolutionary feature and for each 
Z - the impact of the EOS change. A comparison of these values 
with the one of the "Impact" column in Table [2] shows that the 
EOS change accounts for about 10% to 20% of the total impact 
of all the other physical input. The only exception is the helium 
core mass: in this case the impact of the EOS variation accounts 
for 20% to 40% of the total impact of the other physical input. 



Table 12. Uncertainty on the age inferred from the BTO luminosity. 
The uncertainty is evaluated around 12 Gyr. 



z 


Uncertainty (Gyr) 


0.0001 


0.325 


0.0005 


0.345 


0.0010 


0.355 


0.0060 


0.345 


0.0100 


0.365 


0.0200 


0.415 



Article number, page 6 ofllOl 



Valle, G. et al.: Cumulative physical uncertainty in modern stellar models (RN) 



Table A.l. Impact of equation of state change from OPAL to FreeEOS 
for the examined evolutionary features 



z 


log Lbto 


?H 


log Lap 


M c He 


log L HB 




(dex) 


(Gyr) 


(dex) 


(M ) 


(dex) 


0.0001 


0.0017 


0.053 


0.0022 


-0.00095 


-0.0039 




(6%) 


(13%) 


(7%) 


(17%) 


(8%) 


0.0005 


0.0025 


0.049 


0.0005 


-0.00117 


-0.0037 




(10%) 


(11%) 


(2%) 


(22%) 


(10%) 


0.0010 


0.0031 


0.057 


-0.0006 


-0.00132 


-0.0046 




(13%) 


(12%) 


(2%) 


(28%) 


(12%) 


0.0060 


0.0028 


0.096 


-0.0020 


-0.00149 


-0.0070 




(14%) 


(13%) 


(7%) 


(36%) 


(17%) 


0.0100 


0.0023 


0.110 


-0.0027 


-0.00142 


-0.0030 




(12%) 


(13%) 


(9%) 


(38%) 


(7%) 


0.0200 


0.0031 


0.192 


-0.0027 


-0.00115 


-0.0013 




(18%) 


(18%) 


(9%) 


(33%) 


(3%) 



Notes. In parenthesis: percentage influence of EOS change with respect 
to the values in the column "Impact" of Table|2] 
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Table A.3. Estimated coefficients for the linear model in Eq.[T]of BTO 
luminosity (dex), for different Z values. 



A 

A (pp) 
A ( 14 N) 

A (h) 

Pa (vg) 



A 

Pi (PP) 

A ( 14 N) 

Pi (h) 

Pa (vg) 



A) 
Ai (pp) 

A2 ( 14 N) 

A (fe) 
At (v d ) 



A) 

Ai (pp) 
Pi ( 14 N) 

A (*r) 

A (v d ) 



A 

Ai (pp) 
Pi ( 14 N) 

A (*r) 

A (vd) 



A 

A (pp) 
A ( 14 N) 
A (fe) 
A (vd) 



Estimate 
(dex) 



Std. Error 
(dex) 



Impact 
(dex) 



Importance 



Z = 0.0001 



1.1455 

0.0581 

-0.0338 

-0.4050 

-0.0160 



0.0026 
0.0020 
0.0006 
0.0012 
0.0004 



0.0017 
-0.0034 
-0.0202 
-0.0024 



Z = 0.0005 



0.9593 

0.0730 

-0.0312 

-0.3411 

-0.0124 



0.0027 
0.0021 
0.0006 
0.0012 
0.0004 



0.0022 
-0.0031 
-0.0171 
-0.0019 



Z = 0.001 



0.8775 

0.0735 

-0.0315 

-0.3229 

-0.0124 



0.0026 
0.0020 
0.0006 
0.0012 
0.0004 



0.0022 
-0.0032 
-0.0161 
-0.0019 



Z = 0.006 



0.6031 

0.0650 

-0.0278 

-0.2760 

-0.0098 



0.0028 
0.0021 
0.0006 
0.0012 
0.0004 



0.0020 
-0.0028 
-0.0138 
-0.0015 



Z = 0.01 



0.5144 

0.0620 

-0.0250 

-0.2737 

-0.0074 



0.0030 
0.0023 
0.0007 
0.0014 
0.0005 



0.0019 
-0.0025 
-0.0137 
-0.0011 



Z = 0.02 



0.3715 

0.0571 

-0.0248 

-0.2484 

-0.0044 



0.0025 
0.0019 
0.0006 
0.0011 
0.0004 



0.0017 
-0.0025 
-0.0124 
-0.0007 



0.086 
0.167 
1.000 
0.119 



0.128 
0.183 
1.000 
0.109 



0.137 
0.195 
1.000 
0.115 



0.141 
0.202 
1.000 
0.106 



0.136 
0.183 
1.000 
0.081 



0.138 
0.200 
1.000 
0.053 



Notes. In the first two columns: least-squares estimates of the regres- 
sion coefficients and their errors; third column: physical impact of the 
variation Ap, of the various inputs; last column: relative importance of 
the impact with respect to the most important one. 



Table A.4. Estimated coefficients for the linear model in Eq.[T]of the 
central hydrogen exhaustion time (Gyr), for different Z values. The 
column legend is the same as in Table lA.3l 



A 
A (pp) 

A ( 14 N) 

A (*r) 

Pa (va) 



A 

A (pp) 
A ( 14 N) 

A (kr) 

Pa (v d ) 



A 
A (pp) 

A ( 14 N) 

A (kr) 

Pa (v d ) 



A 

A (pp) 
A ( 14 N) 

A (*r) 

A (v d ) 



A 

A (pp) 
A ( 14 N) 
A (h) 

Pa (v d ) 



A 

A (pp) 
A ( 14 N) 

A (*r) 

A (vd) 



Estimate 
(Gyr) 



Std. Error 
(Gyr) 



Impact 
(Gyr) 



Importance 



Z = 0.0001 



-0.4596 
0.5368 

-0.0776 
7.4100 

-0.1982 



0.0115 
0.0089 
0.0027 
0.0052 
0.0017 



0.016 
-0.008 

0.370 
-0.030 



Z = 0.0005 



-0.7055 
0.6674 

-0.0965 
7.7301 

-0.2115 



0.0105 
0.0081 
0.0024 
0.0047 
0.0016 



0.020 
-0.010 

0.387 
-0.032 



Z = 0.001 



-0.8548 
0.8001 

-0.1047 
8.1423 

-0.2379 



0.0130 
0.0101 
0.0030 
0.0059 
0.0020 



0.024 
-0.011 

0.407 
-0.036 



Z = 0.006 



-2.2543 
1.7872 
-0.1437 
11.5876 
-0.4341 



0.0306 
0.0237 
0.0071 
0.0139 
0.0046 



0.054 
-0.014 

0.579 
-0.065 



Z = 0.01 



-2.4268 
2.2898 
-0.1721 
13.2384 
-0.6116 



0.0421 
0.0325 
0.0097 
0.0191 
0.0064 



0.069 
-0.017 

0.662 
-0.092 



Z = 0.02 



-3.1791 
3.2350 
-0.2074 
16.3678 
-0.9817 



0.0600 
0.0464 
0.0139 
0.0272 
0.0091 



0.097 

-0.021 

0.818 

-0.147 



0.043 
0.021 
1.000 
0.080 



0.052 
0.025 
1.000 
0.082 



0.059 
0.026 
1.000 
0.088 



0.092 
0.025 
1.000 
0.112 



0.104 
0.026 
1.000 
0.139 



0.119 
0.025 
1.000 
0.180 
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Table A.5. Estimated coefficients for the linear model in Eq.[2]of RGB 
tip luminosity (dex), for different Z values. The column legend is the 
same as in TablelA.3l 



Table A.6. Estimated coefficients for the linear model in Eq.[2]of helium 
core mass Mf e (M Q ) for different Z values. The column legend is the 
same as in TablelA.3l 





Estimate 


Std. Error 


Impact 


Importance 




Estimate 


Std. Error 


Impact 


Importance 




(dex) 


(dex) 


(dex) 






(Mo) 


(M e ) 


(Mo) 








Z = 0.0001 










Z = 0.0001 






ft 


3.3628 


0.0120 






ft 


0.4839 


0.0042 






ft (pp) 


-0.0086 


0.0077 


-0.0003 


0.028 


ft (PP) 


-0.0013 


0.0027 


-0.0000 


0.023 


ft ( 14 N) 


0.0377 


0.0023 


0.0038 


0.405 


ft ( 14 N) 


-0.0061 


0.0008 


-0.0006 


0.344 


ft (Ar) 


-0.1858 


0.0046 


-0.0093 


1.000 


ft (Ar) 


0.0147 


0.0016 


0.0007 


0.411 


Pa (vd) 


0.0041 


0.0015 


0.0006 


0.066 


ft (Vd) 


0.0022 


0.0005 


0.0003 


0.185 


ft (3a) 


-0.0402 


0.0011 


-0.0080 


0.866 


ft (3a) 


-0.0089 


0.0004 


-0.0018 


1.000 


ft(v) 


0.1552 


0.0057 


0.0062 


0.668 


ft(v) 


0.0323 


0.0020 


0.0013 


0.724 


Pi (*b) 


-0.0676 


0.0047 
Z = 0.0005 


-0.0034 


0.364 


ft (Ac) 


-0.0143 


0.0016 
Z = 0.0005 


-0.0007 


0.402 


Po 


3.5049 


0.0109 






ft 


0.4806 


0.0021 






Pi (pp) 


-0.0062 


0.0070 


-0.0002 


0.016 


ft (PP) 


-0.0037 


0.0013 


-0.0001 


0.069 


Pi ( 14 N) 


0.0271 


0.0021 


0.0027 


0.231 


ft ( 14 N) 


-0.0057 


0.0004 


-0.0006 


0.353 


Pi (Ar) 


-0.2349 


0.0041 


-0.0117 


1.000 


ft (Ar) 


0.0186 


0.0008 


0.0009 


0.573 


Pa (v d ) 


0.0037 


0.0014 


0.0006 


0.047 


ft (v d ) 


0.0021 


0.0003 


0.0003 


0.196 


ft (3a) 


-0.0362 


0.0010 


-0.0072 


0.616 


ft (3a) 


-0.0081 


0.0002 


-0.0016 


1.000 


ft(v) 


0.1258 


0.0051 


0.0050 


0.428 


ft(v) 


0.0242 


0.0010 


0.0010 


0.596 


A (Arc) 


-0.0693 


0.0042 
Z = 0.001 


-0.0035 


0.295 


ft (Ac) 


-0.0144 


0.0008 
Z = 0.001 


-0.0007 


0.442 


Po 


3.5569 


0.0101 






ft 


0.4749 


0.0017 






Pi (PP) 


-0.0069 


0.0065 


-0.0002 


0.017 


ft (PP) 


-0.0017 


0.0011 


-0.0001 


0.033 


A ( 14 N) 


0.0256 


0.0019 


0.0026 


0.206 


ft ( 14 N) 


-0.0048 


0.0003 


-0.0005 


0.312 


Pi (h) 


-0.2478 


0.0039 


-0.0124 


1.000 


ft (Ar) 


0.0172 


0.0006 


0.0009 


0.555 


Pa (Vd) 


0.0039 


0.0013 


0.0006 


0.048 


ft (v d ) 


0.0020 


0.0002 


0.0003 


0.194 


ft (3a) 


-0.0351 


0.0010 


-0.0070 


0.566 


ft (3a) 


-0.0077 


0.0002 


-0.0015 


1.000 


ft(v) 


0.1171 


0.0048 


0.0047 


0.378 


ft(v) 


0.0230 


0.0008 


0.0009 


0.596 


ft (Ac) 


-0.0719 


0.0039 
Z = 0.006 


-0.0036 


0.290 


ft (Ac) 


-0.0120 


0.0007 
Z = 0.006 


-0.0006 


0.388 


ft 


3.6840 


0.0024 






ft 


0.4779 


0.0004 






ft (pp) 


-0.0020 


0.0015 


-0.0001 


0.004 


ft (PP) 


-0.0013 


0.0003 


-0.0000 


0.027 


ft ( 14 N) 


0.0277 


0.0005 


0.0028 


0.185 


ft ( 14 N) 


-0.0049 


0.0001 


-0.0005 


0.338 


A (Ar) 


-0.2997 


0.0009 


-0.0150 


1.000 


ft(Ar) 


0.0131 


0.0001 


0.0007 


0.456 


/?4 (Vd) 


0.0021 


0.0003 


0.0003 


0.021 


ft (v d ) 


0.0020 


0.0001 


0.0003 


0.209 


A (3a) 


-0.0303 


0.0002 


-0.0061 


0.404 


ft (3a) 


-0.0072 


0.0000 


-0.0014 


1.000 


ft(v) 


0.0816 


0.0011 


0.0033 


0.218 


ft(v) 


0.0158 


0.0002 


0.0006 


0.440 


Pi (*b) 


-0.0574 


0.0009 
Z=0.01 


-0.0029 


0.192 


ft (Ac) 


-0.0119 


0.0002 
Z = 0.01 


-0.0006 


0.412 


ft 


3.7028 


0.0132 






ft 


0.4782 


0.0015 






A (pp) 


-0.0048 


0.0085 


-0.0001 


0.010 


ft (PP) 


-0.0012 


0.0010 


-0.0000 


0.026 


ft ( 14 N) 


0.0267 


0.0025 


0.0027 


0.188 


ft ( 14 N) 


-0.0046 


0.0003 


-0.0005 


0.331 


ft (A r ) 


-0.2845 


0.0050 


-0.0142 


1.000 


ft (Ar) 


0.0095 


0.0006 


0.0005 


0.340 


ft (Vd) 


0.0054 


0.0017 


0.0008 


0.057 


ft (v d ) 


0.0018 


0.0002 


0.0003 


0.195 


ft (3a) 


-0.0296 


0.0013 


-0.0059 


0.416 


ft (3a) 


-0.0070 


0.0001 


-0.0014 


1.000 


ft(v) 


0.0641 


0.0062 


0.0026 


0.180 


ft(v) 


0.0145 


0.0007 


0.0006 


0.413 


ft (Ac) 


-0.0616 


0.0051 
Z = 0.02 


-0.0031 


0.216 


ft (Ac) 


-0.0104 


0.0006 
Z = 0.02 


-0.0005 


0.372 


ft 


3.7215 


0.0116 






ft 


0.4769 


0.0015 






ft (PP) 


0.0000 


0.0075 


0.0000 


0.000 


ft (PP) 


-0.0015 


0.0010 


-0.0000 


0.031 


ft ( 14 N) 


0.0249 


0.0022 


0.0025 


0.167 


ft ( 14 N) 


-0.0041 


0.0003 


-0.0004 


0.292 


ft (Ar) 


-0.2988 


0.0044 


-0.0149 


1.000 


ft (Ar) 


0.0070 


0.0006 


0.0003 


0.249 


ft (Vd) 


0.0030 


0.0015 


0.0004 


0.030 


ft (v d ) 


0.0023 


0.0002 


0.0003 


0.244 


ft (3a) 


-0.0280 


0.0011 


-0.0056 


0.374 


ft (3a) 


-0.0070 


0.0001 


-0.0014 


1.000 


ft(v) 


0.0642 


0.0055 


0.0026 


0.172 


ft(v) 


0.0117 


0.0007 


0.0005 


0.334 


ft (Ac) 


-0.0553 


0.0045 


-0.0028 


0.185 


ft (Ac) 


-0.0103 


0.0006 


-0.0005 


0.369 
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Table A.7. Estimated coefficients for the linear model in Eq. [2] of the 
ZAHB luminosity log Lhb at log T e g = 3.83 (dex) for different Z values. 
The column legend is the same as in Table lA.3l 





Estimate 


Std. Error 


Impact 


Importance 




(dex) 


(dex) 


(dex) 




Z = 0.0001 


A 


2.3509 


0.0168 






A (pp) 


-0.0150 


0.0108 


-0.0004 


0.017 


A ( 14 N) 


-0.0594 


0.0032 


-0.0059 


0.227 


A (k) 


-0.5224 


0.0064 


-0.0261 


1.000 


Pa (v d ) 


-0.0113 


0.0021 


-0.0017 


0.065 


A (3a) 


-0.0237 


0.0016 


-0.0047 


0.182 


AM 


0.1267 


0.0079 


0.0051 


0.194 


Pi (k c ) 


-0.0539 


0.0065 


-0.0027 


0.103 


Z = 0.0005 


A 


2.1457 


0.0091 






Pi (pp) 


-0.0019 


0.0058 


-0.0001 


0.003 


Pi ( 14 N) 


-0.0015 


0.0017 


-0.0002 


0.007 


Aft) 


-0.4209 


0.0035 


-0.0210 


1.000 


Pa (vd) 


-0.0093 


0.0012 


-0.0014 


0.066 


A (3a) 


-0.0399 


0.0009 


-0.0080 


0.380 


AM 


0.0712 


0.0043 


0.0028 


0.135 


Pi (k c ) 


-0.0420 


0.0035 


-0.0021 


0.100 


Z = 0.001 


A 


2.0878 


0.0121 






Ai (pp) 


0.0163 


0.0078 


0.0005 


0.024 


A ( 14 N) 


0.0106 


0.0023 


0.0011 


0.052 


A (*r) 


-0.4126 


0.0046 


-0.0206 


1.000 


A (v d ) 


-0.0078 


0.0015 


-0.0012 


0.057 


A (3a) 


-0.0446 


0.0012 


-0.0089 


0.432 


AM 


0.0623 


0.0057 


0.0025 


0.121 


A ft) 


-0.0415 


0.0047 


-0.0021 


0.101 


Z = 0.006 


A 


2.0339 


0.0045 






A (pp) 


0.0040 


0.0029 


0.0001 


0.006 


A ( 14 N) 


0.0133 


0.0009 


0.0013 


0.062 


A ft) 


-0.4304 


0.0017 


-0.0215 


1.000 


A (v d ) 


-0.0109 


0.0006 


-0.0016 


0.076 


A (3a) 


-0.0579 


0.0004 


-0.0116 


0.538 


AM 


0.0556 


0.0021 


0.0022 


0.103 


A ft) 


-0.0427 


0.0017 


-0.0021 


0.099 


Z=0.01 


A 


1.9961 


0.0067 






A. (pp) 


0.0045 


0.0043 


0.0001 


0.006 


A ( 14 N) 


0.0362 


0.0013 


0.0036 


0.161 


A ft) 


-0.4489 


0.0025 


-0.0224 


1.000 


A (vd) 


-0.0168 


0.0009 


-0.0025 


0.112 


A (3a) 


-0.0583 


0.0006 


-0.0117 


0.520 


AM 


0.0309 


0.0032 


0.0012 


0.055 


A ft) 


-0.0215 


0.0026 


-0.0011 


0.048 


Z = 0.02 


A 


1.9478 


0.0056 






A. (pp) 


0.0084 


0.0036 


0.0003 


0.011 


A ( 14 N) 


0.0425 


0.0011 


0.0042 


0.185 


A ft) 


-0.4588 


0.0021 


-0.0229 


1.000 


A (vd) 


-0.0189 


0.0007 


-0.0028 


0.124 


A (3a) 


-0.0603 


0.0005 


-0.0121 


0.526 


AM 


0.0211 


0.0027 


0.0008 


0.037 


A ft) 


-0.0172 


0.0022 


-0.0009 


0.037 



Table A.8. Stellar models computed for isochrones construction at dif- 
ferent Z values. 

~~Z Stellar models (M G ) 



0.0001 



0.0005 



0.0010 



0.0060 



0.0100 



0.0200 



0.50, 0. 
0.82, 0. 
0.50, 0. 
0.82, 0. 
0.50, 0. 
0.85,0. 
0.50, 0. 
0.92, 0. 
0.70, 0. 
0.97, 1. 
0.70, 0. 
0.97, 1. 



60, 0.70, 0.72, 
85, 0.87, 0.90, 
60, 0.70, 0.72, 
85, 0.87, 0.90, 
60, 0.70, 0.75, 
87, 0.90, 0.92, 
60, 0.70, 0.80, 
95, 1.00, 1.05, 
80, 0.85, 0.87, 
00, 1.02, 1.05, 
80, 0.85, 0.87, 
00, 1.02, 1.05, 



0.75, 0.77, 0.80, 

0.95 

0.75, 0.77, 0.80, 

0.95 

0.77, 0.80, 0.82, 

0.95 

0.85, 0.87, 0.90, 

1.10 

0.90, 0.92, 0.95, 

1.10 

0.90. 0.92, 0.95, 

1.10 



